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Abstract 

Backgrounds: The process of amyloid proteins aggregation causes several human neuropathologies. In some 
cases, e.g. fibrillar deposits of insulin, the problems are generated in the processes of production and purification of 
protein and in the pump devices or injectable preparations for diabetics. Experimental kinetics and adequate 
modelling of chemical inhibition from amyloid aggregation are of practical importance in order to study the viable 
processing, formulation and storage as well as to predict and optimize the best conditions to reduce the effect of 
protein nucleation. 

Results: In this manuscript, experimental data of insulin, AP42 amyloid protein and apomyoglobin fibrillation from 
recent bibliography were selected to evaluate the capability of a bivariate sigmoid equation to model them. The 
mathematical functions (logistic combined with Weibull equation) were used in reparameterized form and the 
effect of inhibitor concentrations on kinetic parameters from logistic equation were perfectly defined and 
explained. The surfaces of data were accurately described by proposed model and the presented analysis 
characterized the inhibitory influence on the protein aggregation by several chemicals. Discrimination between true 
and apparent inhibitors was also confirmed by the bivariate equation. EGCG for insulin (working at pH = 7.4/T = 37°C) 
and taiwaniaflavone for A(342 were the compounds studied that shown the greatest inhibition capacity. 

Conclusions: An accurate, simple and effective model to investigate the inhibition of chemicals on amyloid protein 
aggregation has been developed. The equation could be useful for the clear quantification of inhibitor potential of 
chemicals and rigorous comparison among them. 

Keywords: Protein fibrillation kinetics. Mathematical modeling. Sigmoid bivariate equation 



Background 

The aggregation and fibrillation of proteins has been com- 
monly associated with numerous degenerative disorders 
in humans including Alzheimer's, Parkinson's, prion's, 
diabetes type II and Huntington's diseases [1-5]. These 
proteins, called amylodogenics, are also involved in pro- 
duction, purification and formulation interferences of 
biotechnology and pharmacologic processes [6,7]. The 
phenomenon of protein aggregation, led by nucleation- 
dependent polymerization mechanism, consists in the 
formation of amyloid fibrils formulated with p-sheet 
structures in cross-|3-sheet arrangement producing in- 
soluble aggregates [8-10]. For example, insulin peptide 
also generates fibrillar structures under specific conditions 
such as low pH, high temperatures or organic solvents 
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dilution [11-13]. Those phase transformations are not 
observed in vivo but are very common, for example, in 
the clinical preparations of insulin for diabetics [14]. 
Human amyloid proteins (Ap) are peptides of rather 39- 
42 residues. Ap40 contains 40 amino acids and Ap42 is 
the major isoform in the Ap peptides with 42 residues 
polypeptide chain and it is the responsible of amyloid 
plaques generated in Alzheimer's disorder [15,16]. 

In general, much of the success in the future application 
of chemicals to inhibit the in vivo formation of such fibrils 
is dependent on the correct modelling and evaluation of 
in vitro kinetic data in order to establish protocols of ac- 
tion, effectiveness of molecules and dose-optimization of 
compounds. The kinetic description of amyloid protein 
aggregation based on mechanistic and thermodynamic ap- 
proaches have been extensively studied [17-22] and an 
ample number of plausible mechanisms of nucleation and 
fibrillation have been proposed [23,24]. So remarkable, 
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two new and similar proposal: "Ockham's razor'Vminima- 
listic and Crystallization-like Model, have been recently 
developed. Both have solid biophysical basis and they were 
successfully applied to describe and explain the ex- 
perimental data of different amyloid protein aggrega- 
tion [25-28]. However, the combined effect of anti- 
aggregation protein agents, e.g. surfactants, osmolytes and 
food additives, on time-dependent responses has not been 
modelled by that way or by means of empirical equations. 

Although the use of empirical sigmoid equations, 
mainly the logistic equation, does not provide a direct 
explanation of the molecular steps that underlie in the 
generation of fibrils, it is a robust tool to examine pro- 
tein aggregation kinetic data and to address all the 
phases of the process [4,29,30]. In addition, other au- 
thors assimilated the parameters from logistic [20] to 
those obtained by an autocatalytic mechanism [23,27]. 
Nonetheless, that equation is always formulated with- 
out the parameters (fibrillation rate and lag phase) in 
an explicit form hindering the estimation of their stat- 
istical error. 

In the present work, the capability of fit and experi- 
mental data predictability of a sigmoid bivariate model 
that simulates the growth of aggregation process on dif- 
ferent proteins along with the effects of inhibitory che- 
micals on the kinetic parameters is explored in selected 
cases obtained from the literature. The results reveal its 
efficacy and validity to analyze the most relevant param- 
eters that describe geometrically and macroscopically 
the mentioned process. 

Methods 

Experimental data 

Amyloid protein aggregation data were collected from 
results previously reported in the bibliography and digi- 
tized from the published curves using GetData Graph 
Digitizer 2.24. The kinetics of insulin inhibition induced 
by (-)-epigallocatechin-3-gallate (EGCG) were selected 
from Wang et al. [31], methylglyoxal effects were col- 
lected from Oliveira et al. [32] and those produced by 
l,2-diheptanoyl-sn-glycero-3-phosphocholine (di-C7-PC) 
were described in Wang et al. [33]. Two conditions of 
EGCG affecting to fibrillation kinetics were used: EGCG_1 
(case 1) studied at pH = 2.0/T = 60°C and EGCG_2 (case 
2) at pH = 7.4/T = 37°C. 

On the other hand, the aggregation kinetics of Ap42 
amyloid protein inhibited by apigenin and taiwaniafla- 
vone were selected from Thapa et al. [16] and the data af- 
fected by ectoine and hydroxyectoine from Kanapathipillai 
et al. [15]. Finally, apomyoglobin fibrillation expe- 
riments were published in Vilasi et al. [34]. The data- 
sets were obtained by the two most common methods 
used to probe amyloid formation in vitro, the incre- 
ment in light scattering of the protein solution due to 



insolubilization, and the increase in ThT fluorescence due 
to amyloid binding. 

Mathematical modelling 

The model developed to simulate the process of aggrega- 
tion and hence insulin fibrillation was defined by a bivari- 
ate equation. Such model is based on the combination of 
WeibuU function as chemical-concentration model [35,36] 
modifying the most important parameters of the repara- 
meterized logistic equation [37] used for aggregation de- 
scription. This expression (1) has been successfully used, 
in recent works, to evaluate the inhibitory effect of organic 
acids and heavy metals on the growth of various bacteria 
[38,39]. Its mathematical form is as follows (see also 
Appendix section): 



X 



1 + exp 2 



; where : 



(1) 



A. = A<^ 1 + Kx 



l-exp(-l„2(-£)* 

1- exp [ - ln2 

where, is the maximum aggregation rate, is the 
maximum aggregation growth, A is the lag phase and C 
is the chemical concentration. The meanings of other 
symbolic notations as well as the corresponding units 
are summarized in Table 1. In the experimental data 
chosen, the dependent variable of response or time- 
dependent signals (X) to detect amyloid protein aggrega- 
tion were absorbance at 600 nm (data obtained from 
Wang et al. [31]), relative ThT fluorescence intensity (%) 
(data from Wang et al. [33]) and ThT fluorescence inten- 
sity at 482 nm or 490 nm (data obtained from Thapa 
et al. [16]; Kanapathipillaia et al. [15]; Oliveira et al. [32] 
and Vilasi et al. [34]). 

Additionally, a global parameter {ECso.t) was also se- 
lected for the overall description of chemical effects on 
fibrillation growth studies according the algebraic steps 
previously described [38]. This parameter was defined as 
the chemical concentration (in mM) that reduces the ag- 
gregation by 50% compared to that produced by the 
control without agent at time (t) which also reduces the 
maximum aggregation by 50%. 

Numerical methods and statistical analysis 

The fitting procedures and parametric estimates from 
the experimental results were performed by minimizing 
the sum of quadratic differences between the observed 
and model-predicted values using the nonlinear least- 
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Amyloid aggregation growtli measured as absorbance at 500 nm, relative ThT fluorescence intensity 
{%) and ThT fluorescence intensity at 482 nm or 490 nm. Units: absorbance units (AU) or (%). 

Time. Units: h or d 

Maximum aggregation growth. Units: AU or % 
Maximum aggregation rate. Units: AU h"', AU d"^ or % h~' 
Lag phase. Units: h or d 

Maximum insulin aggregation affected by chemical agent. Units: AU or % 
Maximum insulin aggregation rate affected by chemical agent Units: AU h"\ AU or % h~^ 
Lag phase affected by chemical agent. Units: h or d 

Concentration of chemical agent. Units: mM or ^iM 
Maximum response affecting on X^. Dimensionless 

Concentration corresponding to the semi-maximum response affecting on X^. Units: mM or 
Shape parameter affecting on X^. Dimensionless 
Maximum response affecting on v^. Dimensionless 

Concentration corresponding to the semi-maximum response afl'ecting on v^. Units: mM or pM 
Shape parameter affecting on Vm. Dimensionless 
Maximum response affecting on X. Dimensionless 

Concentration corresponding to the semi-maximum response affecting on A. Units: mM or pM 
Shape parameter affecting on A. Dimensionless 



Table 1 Symbolic notations used and corresponding units 

nsulin aggregation kinetics measured by 
absorbance or fluorescence 

X: 
t : 

'■ 

Vm ■■ 

A : 
Vm- ■■ 

Concentration effects on insulin aggregation kinetics 
C: 
Kx- 

: 

Ox : 
: 

: 

Qv, : 
K,: 

: 

: 



squares (quasi-Newton) method provided by the 'Solver 
macro from Microsoft Excel spreadsheet. The confi- 
dence intervals of the best-fit values for the parametric 
estimates (Student's t test, a = 0.05), consistency of the 
mathematical models (Fisher's F test; p < 0.05) and covari- 
ance and correlation matrices were calculated using the 
'SolverAid' macro, which is freely available from Levie's 
Excellaneous website http://www.bowdoin.edu/~rdelevie/ 
exellaneous/. An example of the type of Excel spreadsheet 
used for modeling (Ab42amyloid apigenin case) was pro- 
vided as Additional file 1. These statistical procedures and 
residual analysis (Durbin- Watson test) were confirmed 
and evaluated, respectively, by DataFit 9 (Oakdale Engi- 
neering, Oakdale, PA). Moreover, bias and accuracy fac- 
tors of fittings (Bf and Ajj respectively) were calculated 
according to the expressions reported by Ross [40] . 

Results and discussion 

Characteristics and simulations of bivariate model 

In the description of amyloidogenic fibrillation growths, 
the logistic equation used to formalize their kinetic pro- 
files is always formulated through an explicit expression 
based on the time required to reach 50% of the maximal 
aggregation {xq) and a time parameter (t), although the 
initial lag phase and the aggregation rate are commonly the 
most important parameters that need to be calculated [30]. 
In this work, the logistic equation was reparameterized to 



make the lag phase (A) and the maximum aggregation rate 
{vm)> which represents the overall aggregation process rate, 
explicit (Figure 1). This last coefficient more adequately 
characterizes the reaction rate of fibrillation than the par- 
ameter conventionally used {kapp, defined as specific or ap- 
parent rate) because it is less sensitive to experimental 
error and because it is a more efficient to describe the glo- 
bal rate of the kinetic process [41]. Additionally, the cal- 
culated from the experimental data obtained by Sabate 
et al. [42] for yeast prion proteins follows the Arrhenius de- 
pendence with temperature (data not shown). 

In contrast, the fit using the reparameterized functions 
can be used to easily calculate the confidence intervals 
of the parameters. The algebraic steps required to obtain 
the corresponding reparameterization of A and v„ from 
the logistic equation are detailed in the Appendix. 
Thus, the parameters from this equation define all of the 
aggregation kinetic phases [19]: pre-nucleation, which is 
characterized by A, post-nucleation, which is represented 
by or x, and elongation, which is determined by X^,. 

The absolute correlation between the lag phase and the 
aggregation growth was recently demonstrated in several 
sets of protein data [43]. Indeed, this relationship is obvi- 
ous because the specific or apparent rate (k^pp or kg, de- 
pending on the authors) is inversely proportional to the 
lag phase based the following expression when the logistic 
model (1) is applied [44]: 
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Figure 1 Left, Graphical description of the kinetic parameters (X„, X, v„ and t) from the logistic equation (A.8) and the corresponding 
aggregation phases: pre-nucleation, post-nucleatlon and elongation. Right, Simulations of the most common profiles for the parameters 
(Xm; Vm-), affected by chemical concentration, using the Weibull equations (A.15). 



^app 



Figure 2 shows an illustrative set of simulations that 
were developed under the numerical conditions de- 
scribed in Table 2. In these simulations, all of the pos- 
sible sigmoid effects of the chemical concentration on 
the logistic parameters are defined. WeibuU's equation, 
which was adequately configured to simulate the dose- 
response trends, is a mathematical tool that provides ex- 
cellent predictions for varied experimental profiles from 
different scientific fields [45,46]. In the present con- 
text, the flexibility of this equation to describe different 
chemical-concentration relationships as linear, sigmoid, 
exponential, or hyperbolic curves depending on the nu- 
merical values of the parameters is particularly valuable. 
The combination of both of these mathematical resources 
(1) generates a wide variety of theoretical casuistries 
that will be verified using the experimental data se- 
lected below. 

Inhibitory effect of EGCG on insulin aggregation 

The representation of the experimental data from se- 
lected cases and surfaces predicted by equation (1) are 
depicted in Figure 3A and B. The results of the param- 
eter estimations and the statistical analysis of the fit for 
the EGCG influence on insulin fibrillation are summa- 
rized in Table 3. Both datasets were accurately fitted 
with equation (1) to obtain 7?f values greater than 0.98. 
Moreover, the /^-values obtained from Fisher's F test in- 
dicated the complete consistency and robustness of the 
equation to the adjustment of the experimental patterns. 
The residual analysis through the Durbin- Watson test 
revealed a lack of residual autocorrelation and therefore 
a random distribution (data not shown). 



At an acidic pH and a high temperature, the effect of 
EGCG on the kinetic parameters was not significant for 
the maximum aggregation rate but statistically signifi- 
cant for and the lag phase. As can be observed in 
Figure 3A, the slopes of the sigmoid curves were parallel 
to the influence of the EGCG concentration, and the 
surface can be assimilated to that simulated in Figure 2D 
but with less effect on X,„. These results show a longer 
pre-nucleation phase, which indicates that EGCG blocks 
the formation of the seeding nuclei without changing 
the fibril elongation rate (post-nucleation). This behav- 
iour is consistent with the definition reported by Martins 
[26] who denominates the compounds that modify kin- 
etic parameters (for instance in the present case, and/ 
or X) as true inhibitors. 

In contrast, different surface and parametric responses 
were observed for EGCG_2 up to a concentration of 0.2 
mM polyphenol with only significant modifications in 
the maximum aggregation growth parameter (Figure 3B). 
Under experimental conditions of higher pH and lower 
T, EGCG acts as an apparent inhibitor that only alters 
the thermodynamic properties of aggregation process 
[26]. Higher concentrations of EGCG (from 0.35 to 0.7 
mM) led to an increase in the fibrillation response more 
than the inhibition of the amyloid process [31]. These 
experimental points were deleted in this mathematical 
analysis because biphasic aggregation data with profiles 
of decreasing and increasing values of cannot be 
fitted by the proposed equation or by any other simple 
equation. 

The global description of chemical inhibition was cal- 
culated by means of a single index: ECso.t (Table 3). This 
parameter can be defined as a summary of all of the ef- 
fects on the aggregation kinetics observed at the time re- 
quired to reach the semi-maximum aggregation growth. 
Thus, the obtained value of this parameter can be used 
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Figure 2 Profiles obtained by simulation under the numerical conditions specified in Table 2 using the equation (1) and that 
representing all the theoretical kinetics of protein aggregation process. A: all the parameters from logistic equation {X^, and A) are 
affected by chemical concentration; B: parameters (X^ and v^) are modified by chemical; C: only the parameter is affected by chemical; 
D: parameters and A) are modified by chemical; E: parameters (i/^ and A) are affected by chemical; F: only the parameter is modified by 
chemical; G: only the parameter A is affected by chemical. In all cases, time (r), aggregation response (X) and chemical concentration (Q are 
simulated with arbitrary units. 



for the comparison and the evaluation of the application 
of chemicals to reduce the protein fibrillation. The results 
demonstrated that the most inhibitory conditions for 
EGCG are established at neutral pH and a lower tem- 
perature (case 2), but these conditions require a longer 
time (r) to achieve the parameter ECso.t- 

Inhibitory effect of dl-C7-PC and methylglyoxal on insulin 
aggregation 

Figure 3C displays the surface described by equation (1) 
and data on the insulin kinetics in the presence of a 



surfactant. The parameter values and the statistical char- 
acteristics of the modeling were also determined (Table 3). 
All of the statistical results revealed good agreement and 
accuracy between the observed and the predicted values, 
the consistency of the equation, and a lack of bias during 
the fitting process. 
Regarding the effects, the two kinetic parameters and 
were significantiy affected by di-C7-PC, whereas the 
lag phase was not modified by the surfactant (Student's t 
test, a = 0.05), as determined in the simulated conditions 
defined by Figure 2B. This result is in agreement with data 
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Table 2 Arbitrary numerical values defined for the 
simulations of Figure 1 (A, B, C, D, E, F and G) according 
to the parameters defined in the equation (1) 

Simulation conditions 



Parameters 




A 


B 


c 


D 


E 


F 


G 


Aggregation growth 




1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 




Vm 


0.25 


0.25 


0.25 


0.25 


0.50 


0.20 


050 




A 


3.00 


3.00 


3.00 


3.00 


4.00 


5.00 


4.00 


Effect on Xm 




1.00 


1.00 


1.00 


1.00 












5.00 


5.00 


5.00 


5.00 










Ox 


2.00 


2.00 


2.00 


2.00 








Effect on Vm 


Kv 


1.00 


1.00 






0.60 


0.60 








4.00 


4.00 






3.00 


8.00 








2.00 


2.00 






2.00 


2.00 




Effect on A 


Ka 


1.00 






1.00 


4.00 




4.00 






2.00 






2.00 


10.00 




10.00 






2.00 






2.00 


2.00 




2.00 



reported for different systems of effector/protein, such as 
copper/p-amyloid [47] or glutathione/hen egg-white lyso- 
zyme, although these data were obtained using a non- 
mathematical approach [48] . The value of ECso.t was simi- 
lar to that obtained for EGCG (case 1) for a lower initial 
concentration of insulin, but a much longer time was re- 
quired to obtain this parametric value (150.68 h). 

Similar inhibitory responses were also observed with 
methylglyoxal; thus, and were the kinetic parame- 
ters inhibited by this aldehyde form of pyruvic acid 
(Table 3 and Figure 3D). Under an initial insulin concen- 
tration (3 g/L) that was higher than the initial EGCG 
concentration, methylglyoxal effectively reduced the ag- 
gregation growth (0.33 mM) within a shorter time (3.9 
h). Both chemicals (di-C7-PC and methylglyoxal), under 
the experimental conditions assayed, showed similar fea- 
tures of true inhibitors [26]. The application of the pro- 
posed methodology may help the design of optimal 
strategies for the reduction of the fibril deposition of in- 
sulin in commercial preparations. 




EGCG (mM) 



120 




di-C7-PC (mM) 



B 




400 



200 




Methylglyoxal (mM) 



Figure 3 Insulin fibrillation kinetics at different concentrations of EGCG, di-C7-PC and methylglyoxal (points) and fittings to equation 
(1) (surface). X: aggregation data measured by absorbance at 600 nm (AU), ThT fluorescence intensity at 482 nm (AU) or ThT fluorescence 
intensity (%). A: experimental data of EGCG obtained at pH = 2.0 and T = 60°C. B: experimental data of EGCG obtained at pH = 7.4 and T=37°C 
C: di-C7-PC data. D: methylglyoxal data. 
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Table 3 Parametric estimates and confidence intervals (a = 0.05) from the equation (1) applied to the aggregation 
insulin increased data influenced by EGCG_1, EGCG_2, di-C7-PC and methylglyoxal concentrations 



Parameters 




EGCG_1 


EGCG_2 


di-C7-PC 


Methylglyoxal 


Aggregation model 


Xm (AU or %) 


1.12±0.03 


1 .28±0.06 


99.85±4.91 


541 .20±1 1.38 




v„, (AU h"^ or% h"^) 


0.18±0.01 


0.25±0.04 


1.55±0.27 


219.68±20.98 




A(h) 


4.78±0.28 


53.34±0.28 


1 24.62±5.86 


2.92±0.12 


Effect on Xm 


K. 


043±0.03 


0.94±0.02 


0.75±0.05 


1 .00±0.44 




(mM) 


0.68±0.08 


0.01 ±0.00 


0.65±0.12 


1.80±1.25 






1 .65±0.45 


0.96±0.20 


4.04±2.37 


0.99±0.23 


Effect on 


Kv 


NS 


NS 


0.86±0.04 


0.99±0.01 




(mM) 


NS 


NS 


0.47±0.06 


0.42±0.07 






NS 


NS 


2.51±235 


0.74±0.08 


Effect on A 


K, 


0.81 ±0.10 


NS 


NS 


NS 




(mM) 


0.50±0.08 


NS 


NS 


NS 




ax 


0.98+0.24 


NS 


NS 


NS 




EQoj (mM) 


0.44 


0.006 


0,47 


0.33 




T(h} 


7.65 


55.02 


1 50.68 


3.92 




p-value 


<0.001 


<0.001 


<0.001 


<0.001 




Bf 


0.88 


0.80 


1.00 


1.04 




Af 


1.23 


1.40 


1.18 


1.15 






0.992 


0.977 


0.984 


0.974 



Statistical values of adjusted coefficient of multiple determination (R^dj) 3nd p-values from Fisher's F-test (a = 0.05). Bf and Af are the bias and accuracy factor, 
respectively. NS: non-significant. 



Inhibitory effect of mono and biflavonoids on A|342 
amyloid protein 

The dependence of the kinetic parameters of fibril- 
lation on the apigenin and taiwaniaflavone concen- 
trations is represented in Figures 4A and B. Both 
response surfaces were similar to that defined by the 
simulation shown in Figure 2A. The three parameters 
{v^, X^, and A) were clearly and significantly modified 
by these chemical doses (Table 4), which show a longer 
pre-nucleation phase, a lower elongation rate (post-nu- 
cleation), and a lower elongation stage. Such responses 
confirmed their chemical features of true inhibitors, 
being the only case in which both kinetic parameters 
{Vm and A) significantly affected aggregation process 
[26]. All of the statistical tests corroborated the perfect 
agreement between the experimental and the theoret- 
ical data {e.g., Radj > 0.99). The bias and accuracy fac- 
tors {Bf and Aj) were also indicative of the good fitting 
quality obtained with model (1). 

The numerical values of ECso,t confirmed the superior 
antifibrillogenic capacity of the biflavonoid taiwaniafla- 
vone {ECso,T = 0-91 l^M) compared with the monofla- 
vonoid apigenin {ECso.r = 2.45 ^M). However, the time t 
was quite similar in both cases. This affirmation is con- 
sistent with the conclusions reported in the original 
work conducted by Thapa et al. [16]. 



Inhibitory effect of osmolytes on AP42 amyloid protein 

Different types of molecules have been recently sug- 
gested as drug candidates for the treatment of neurode- 
generative disorders caused by the anti-aggregation 
properties of amyloid proteins [49] . One of the most inter- 
esting candidates are osmolytes (compounds that protect 
against the osmotic stress), such as the ectoines that are 
obtained from aerobic heterotrophic bacteria [50]. The re- 
sponse surface for hydroxyectoine (Figure 4C) shows a be- 
havior that is similar to that obtained with the simulation 
shown in Figure 2B. Based on the numerical fittings sum- 
marized in Table 4, it was statistically demonstrated that 
osmolytes inhibited the maximum aggregation rate and 
the maximum fibrillation growth. In contrast, no signifi- 
cant changes in the lag phase were observed. Therefore, 
both osmolytes effectively blocked the post-nucleation 
and the elongation stages but not the pre-nuclei forma- 
tion. In addition, ectoine was much more effective as an 
antineurotoxic {ECso,t = 0.96 mM) than hydroxyectoine 
(£Cso,^ = 2.20 mM). 

Inhibition of apomyoglobin fibrils formation by trehalose 

The effect of trehalose on the apomyoglobin aggregation 
kinetics was also studied using model (1), and excellent 
statistical results were obtained from the modeling and 
the description of the experimental data (Figure 4D and 
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Figure 4 Amyloid protein aggregation {nineties at different concentrations of apigenin, ectoine, taiwaniaflavone and trehalose (points) 

and fittings to equation (1) (surface). X: aggregation data measured by ThT fluorescence intensity at 482 nm or 490 nm (AU). A: experimental 

data of apigenin. B: taiwaniaflavone data. C: experimental data of hydroxyectoine. D: data of trehalose affecting apomyoglobin kinetics. 
^ J 



Table 4). As in the previous case, the disaccharide sig- 
nificantly inhibited the two kinetic parameters Vi„ and 
without affecting the pre-nucleation phase (A). 
Nevertheless, its relative efficacy was obtained at a 
much higher concentration (21.43 mM), and much 
longer times (6.5 d) were required than compared with 
the other cases studied. 

Conclusions 

In summary, a general bivariate model that combines 
the logistic equation for the description of kinetics and 
the WeibuU equation for the chemical-concentration 
response has been proposed for the characterization of 
the inhibitory effects produced by several chemicals on 
the growth of the aggregation of amyloid proteins. In 
all cases, the inhibitory effects on the kinetic parame- 
ters were established, and the theoretical response sur- 
faces were in perfect agreement with the selected data. 
In addition, the recent definition of true and apparent 



inhibitors reported by Martins [24] can be also evalu- 
ated and validated. EGCG at pH 7.4 and a temperature 
of 37°C was the best option for the reduction of the 
formation of insulin fibrils {ECso,t = ^ l^M) but the 
shortest time (t) was obtained with methylglyoxal. For 
the Ap42 amyloid protein, the biflavonoid taiwaniafla- 
vone produced the highest inhibition response at the 
lowest dose (0.91 |iM) and the shortest time (0.65 h). 
The bivariate equation was validated using data obtained 
by two methods based on different chemical phenomena. 
Although this model can not define the mechanisms of ac- 
tion of chemical inhibitors on protein aggregation, it pro- 
vides a consistent tool for the comparison of the ability of 
such compounds in the inhibition of the protein fibrilla- 
tion process, regardless of the method used for its deter- 
mination, and is a first step for the optimization of the 
in vitro application of them. Further experiments and cor- 
responding modeling should be done to establish its valid- 
ity for in vivo applications of anti-aggregation chemicals. 
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Table 4 Parametric estimates and confidence intervals (a = 0.05) from the equation (1) applied to the aggregation of 
AP42-amyloid protein increased data influenced by apigenin, ectoine, hidroxyectoine and taiwaniaflavone 
concentrations 



Parameters 




Apigenin 


Ectoine 


IHidroxy ectoine 


Taiwaniaflavone 


Trehalose 


Aggregation model 


X„ (AU) 


1 7.43±0.50 


(5.07±0.23)x 10"^ 


(4.96±0.24) X 1 0"^ 


1555±055 


1 0.49±0.28 




(AU h"' or AU d"') 


16.10±1.49 


(0.24±0.05) X 1 0"^ 


(0.30±0.08) X 1 0"^ 


17.01 ±2.09 


331 ±0.33 




X (h ord) 


0.25±0.03 


1 1 .60±2.26 


13.50±2.31 


0.27±0.06 


5.76±0.14 


Effect on Xm 


Kx 


0.41 ±0.06 


1 .00±0.99 


1 .00±0.90 


0.78±0.05 


0.44±0.03 




(mM or [iM] 


6.52±1.93 


11.82± 11.03 


66.70±65.78 


2.32±0.38 


5158±7.27 




Ox 


0.91 ±0.25 


0.19±0.13 


0.43±0.43 


0.98±0.16 


1.36±0.31 


Effect on Vm 




0.38±0.12 


0.83±0.06 


0.86±0.16 


1.14±0.60 


0.53±0.08 




(mM or (iM) 


9.29±S.35 


0.33±0.32 


1.35 (NS) 


1.88±1.73 


33.22±17.26 






1.35±1.13 


0.36±0.14 


0.34±0.20 


0.40±0.21 


0.66±0.33 


Effect on A 


Ka 


1 .76±0.70 


NS 


NS 


1.27±1.13 


NS 




mx (mM or ^iM) 


2.16±0.29 


NS 


NS 


253±1.45 


NS 




ax 


2.39±0.99 


NS 


NS 


2.33±2.31 


NS 




ECso,T (mM or pM) 


2.45 


0.96 


2.20 


0.91 


21.43 




J (h or d) 


0.75 


1731 


20.52 


0.65 


6.54 




p-value 


<0.001 


<0.001 


<0.001 


<0.001 


<0.001 




Bf 


0.99 


1.04 


1.00 


1.02 


1.03 






1.06 


1.09 


1.18 


1.11 


1.08 






0.995 


0.988 


0.984 


0.990 


0.995 



The case of apomyoglobin affected by trehalose is also shown. Statistical values of adjusted coefficient of multiple determination (f?^^/) arid p-values from Fisher's 
F-test (a = 0.05). S^and /A/- are the bias and accuracy factor, respectively. NS: non-significant. 



Appendix 

Reparameterization of logistic equation 

It is well-known the autocatalytic origin of the logistic 
equation based on the following differential equation: 



dX _ ( X 



(A.l) 



which, integrated between Xq— >X and Q^t gives the ex- 
plicit form of aggregation growth as a function of the 
time: 



t- — T — ^ Inl^™"^" 



The value of aggregation when t = ti is: 



(A.4) 



(A.5) 



The slope in the inflection point (v^) is defined by the 
following operation: 



X 



x„ 



+ exp[ln(^)-/c«ppf] 



(A.2) 



The parameter t is defined as the time required to obtain 
the semimaximum fibrillation growth (when X = X^/2): 



X„kapp exp^ln 



(A.6) 



x„ 



1 + 



1 , fXm-Xo 
(A.3) 



The lag phase (A) or pre-nucleation phase is defined as 
the intersection of the tangent at the inflection point with 
the abscissa (Figure 1): 



The inflection point {t = ti) can be obtained when the 
second derivative from (A.l) is equal to zero and the 
abscissa is isolated [31]: 



R = X{ti) + v„^{L-ti) 
with L=X when 7? = 0: 
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A = ti 



2v„ 



(A.7) 



^app 



Reorganizing terms, two reparameterized functions can 
be defined: 



X 



X 



X„ 



1+ exp|2 + |-(A-i; 



Xn 



(A.8) 



(A.9) 



When parameters from both equations are influenced 
by chemical agent concentration, they can be rewritten 
as follows: 



X 



X 



exp - 



4v„. 



{X.-t) 



x„ 



exp 



(A.10) 



(A.11) 



On the other hand, the calculation of the time t for 
ECso.i is obtained by the following formula: 



X„ 



Xf, 
Xm- ln3 



4 1+ exp[|f^(r.-r) 



4v„ 



(A.12) 



Weibull equation and bivariate model (\) 
The most interesting form of Weibull equation for 
dose-response modelling is expressed as follows [35]: 



Y=Y^\ 1-exp 



ln2 



(A.13) 



where, Y is the response and Y^ the maximum re- 
sponse, m is the concentration corresponding to the 
semi-maximum response, C is the concentration and a 
is a shape parameter related to the maximum slope of 
the response. 

This equation can be modified according with the 
graphical tendencies of chemicals effects (Figure 1) on 
the aggregation parameters [38]: 



d=eo{l-Ye) or e=eo{\ + Ye) 



(A. 14) 



Thus, decrease of and X^ and increase of A and t 
were the characteristic responses observed for the inhib- 
ition of amyloid protein fibrillation kinetics induced by 
chemicals: 



Xm- — X„ < 1 -K. 



r. = ri 1 + Kr 



(A.15) 



When the equations (A.15) are inserting directly on 
equation (A.IO) or (A.ll), the bivariate model (1) is obtained. 

Additional file 



Additional file 1: Excel spreadsheet used for modeling the 
Ab42amyloid apigenin case. 
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